PROGRAM BUILDING C C This program computes the growth in total building energy use and average C building energy intensity in a given region from 2005 to 2100, based on C actions carried out from 2006 to 2100 C C INPUT VARIABLES: C C AREAD(NS,NR,0:NY) Floor area (billion m2) at 5-year intervals, read in from C the data file that is created from the Excel worksheets C EINT(NS1,NF1,NR1) Stock average electricity energy intensity (kWh/m2/yr) C ENEW0(NS,NF,NR) Energy intensity of new buildings in 2005 relative to the C stock average C EREN0(NS,NF,NR) Energy intensity of buildings that were renovated in 2005, C after renovation, relative to the stock average C ENEW1A(NS,NF,NR) Absolute energy intensity (kWh/m2/yr) of new buildings after C stringent standards are achieved C EREN1A(NS,NF,NR) Absolute energy intensity (kWh/m2/yr) of renovated buildings C after stringent standards are achieved C IYR0(NS,NF,NR) Year in which the transition from EINT to ENEW1A and EREN1A C begins (choose IYR0>2005) C IYR1(NS,NF,NR) Year in which the transition from EINT to ENEW1A and EREN1A C is complete (choose IYR1 no greater than 2050) C FRENO(NS,NR) Fraction of initial building stock that will have been C renovated by 2050 C FREPL(NS,NR) Fraction of initial building stock that will have been C replaced by 2050 C C PRIMARY OUTPUT VARIABLES (used in Excel worksheets): C C ETOTAL(NS,NF,NR,NY) Total energy use (TWh/yr) C AREAF(NS,NR,NY) Floor area year by year (billions m2) C C DETAILED OUTPUT VARIABLES (not used in Excel worksheets): C C ARENO(I) Area of initial building stock group from cohort I that is renovated C (and done so in year i) C AREPL(I) Area of initial building stock group from cohort I that is replaced C (and done so in year i) C AREMN(I) Area of initial building stock from cohort i that remains after its C opportunity for renovation or replacement in year i C ANEW(I) Increase in building area in year I C ATYEAR(I) Total building area in year I C EISTRT(I) Initial energy intensity (kWh/m2/yr) of cohort I C ENEW(I) Energy intensity of cohort I after replacement, or of buildings C added to the total floor area year I C EREN(I) Energy intensity of cohort I after renovation C ETYEAR(I) Energy use (TWh/yr) of cohort I after it has been altered or added to C EIYEAR(I) Average energy intensity (kWh/m2/yr) of cohort I after it has been C altered or added to C ETSTRT(I) Energy use (Twh/yr) of cohorts not yet altered or added to C ETOTAL(I) Total energy use (TWh/yr) in year I, summed over all cohorts C C NOTES: C C NS = number of building sectors (NS1=NS+1) C NF = number of energy sources (NF1=NF+1) C NR = number of regions (NR1=NR+1) C NY = number of years C C Sector 1 is the residential sector C Sector 2 is the commercial sector C NS1 is for the total across sectors C C Fuel 1 is electricity C Fuel 2 is fuels C NF1 is for the total across fuels C C The regions are indicated below. NR1 is for the sum over all regions C C Year 0 is 2005. The first year with any change is 2006 C C The existing building stock is divided into a number of groups equal to the C number of years of the first round of renovations. If we assume a distribution C of energy intensities of the existing building stock, the code assumes that C buildings are renovated or replaced starting with the most energy-intensive group, C assuming that to be the oldest group. This assumption is only an approximation. C C Electricity use in residential buildings is largely for consumer goods, and C is assumed to gradually increase in many regions - hence, 2050 is the year C when greater (rather than reduced) energy intensity is achieved in these regions C PARAMETER (NY1=45,NY=95,NR=10,NF=2,NS=2) PARAMETER (NR1=NR+1,NR2=NR+2,NR3=NR+3,NF1=NF+1,NS1=NS+1) REAL EINT(NS1,NF1,NR1) REAL ENEW0(NS,NF,NR),ENEW1(NS,NF,NR),ENEW1A(NS,NF,NR), * EREN0(NS,NF,NR),EREN1(NS,NF,NR),EREN1A(NS,NF,NR) REAL AREAD(NS,NR,0:19),AREAF(NS1,NR3,0:NY), ANEW(NS,NR,NY), * AREPL(NS,NR,NY),ARENO(NS,NR,NY),AREMN(NS,NR,NY),ARED(NS,NR,NY) REAL EISTRT(NS,NF,NR,NY),ETSTRT(NS,NF,NR,NY),ENEW(NS,NF,NR,NY), * EREPL(NS,NF,NR,NY),EREN(NS,NF,NR,NY), * FREPL(NS,NR),FRENO(NS,NR) REAL ATYEAR(NS,NR,NY),ETYEAR(NS,NF,NR,NY),EIYEAR(NS,NF,NR,NY), * ETOTAL(NS1,NF1,NR3,0:NY) C -------!---------!---------!---------!---------!---------!---------!-- INTEGER YEAR(0:NY),IYR0(NS,NF,NR),IYR1(NS,NF,NR) CHARACTER*10 AREAN,UNITS(3),SECTOR(NS1),FUEL(NF),REGION(NR3), * VARIAB(3) DATA AREAN /' Area'/ DATA SECTOR /' Residentl',' Commercl',' Total'/ DATA FUEL /' Electrcty',' Fuels'/ DATA UNITS /' Gm2',' TWh/yr',' kWh/m2/yr'/ DATA VARIAB /' EISTRT',' ENEW',' EREN'/ DATA REGION /'Pacif_OECD','North_Amer',' W_Europe', * ' E_Europe',' Form_SovU','Lat_Americ', * 'Sub-SaharA','MidEt_NAfr','CenPlannAs', *'Other_Asia',' World','OECD+E+FSU','_____Other'/ C Fractions of existing buildings that are replaced or renovated between C 2005 and 2050 DATA ((FREPL(IS,IR),IS=1,NS),IR=1,NR) /NR*0.25,NR*0.25/ DATA ((FRENO(IS,IR),IS=1,NS),IR=1,NR) /NR*0.75,NR*0.75/ C Starting year DATA ISTART /2006/ C Average energy intensities (kWh/m2/yr) of the existing stock in each region C for residential electricity, residential fuels, commercial electricity, c and commercial fuels DATA ((EINT(IS,1,IR),EINT(IS,2,IR),IS=1,NS),IR=1,NR) / C Pacific OECD * 55., 70., 150.,150., C North America * 66.,100., 150.,150., C Western Europe * 35.,140., 100.,130., C Eastern Europe * 20.,180., 100.,200., C Former Soviet Union * 25.,230., 80.,240., C Latin America * 25., 50., 100., 40., C Sub-Saharan Africa * 5.,220., 40., 5., C Middle East and North Africa * 25., 10., 100., 10., C Centrally-planned Asia * 20., 80., 120.,160., C South Asia and Other Pacific Asia * 10., 60., 60., 5./ C ENERGY INTENSITY OF NEW BUILDINGS IN 2005 RELATIVE TO THE STOCK AVERAGE DATA (((ENEW0(IS,IF,IR),IR=1,NR),IF=1,NF),IS=1,NS) / C Electricity and fuels, residential * 5*1.0,5*1.0, 5*0.8,5*1.0, C Electricity and fuels, commercial * 5*1.0,5*1.0, 5*1.0,5*1.0/ C ABSOLUTE ENERGY INTENSITY OF NEW BUILDINGS AFTER ENHANCED STANDARDS (kWh/m2/yr) DATA (((ENEW1A(IS,IF,IR),IR=1,NR),IF=1,NF),IS=1,NS)/ C Electricity and fuels, residential * 5*30.,25.,20.,20.,30.,25., 5*40.,15.,15.,15.,25.,15., C Electricity and fuels, commercial * NR*50., 5*30.,3*10.,20.,10./ C ENERGY INTENSITY OF RENOVATED BUILDINGS IN 2005 RELATIVE TO THE STOCK AVERAGE DATA (((EREN0(IS,IF,IR),IR=1,NR),IF=1,NF),IS=1,NS) / C Electricity and fuels, residential * 5*1.0,5*1.0, 5*1.0,5*1.0, C Electricity and fuels, commercial * 5*1.0,5*1.0, 5*1.0,5*1.0/ C ABSOLUTE ENERGY INTENSITY OF RENOVATED BUILDINGS AFTER ENHANCED STANDARDS (kWh/m2/yr) DATA (((EREN1A(IS,IF,IR),IR=1,NR),IF=1,NF),IS=1,NS)/ C Electricity and fuels, residential * 5*30.,25.,20.,20.,30.,25., 5*50.,25.,25.,15.,30.,25., C Electricity and fuels, commercial * NR*70., 5*40.,3*15.,30.,15./ C Starting and ending years of transition to deep savings DATA (((IYR0(IS,IF,IR),IR=1,NR),IS=1,2),IF=1,2) * /NR*2010,NR*2010,NR*2010,NR*2010/ DATA (IYR1(1,1,IR),IR=1,NR) /3*2030,7*2050/ DATA (IYR1(1,2,IR),IR=1,NR) /NR*2020/ DATA (IYR1(2,1,IR),IR=1,NR) /NR*2020/ DATA (IYR1(2,2,IR),IR=1,NR) /NR*2020/ C Fractional reduction in energy intensity of new or renovated buildings C compared to the year before, during the second round DATA RATE /0.005/ C Convert absolute energy intensities with enhanced standards to C energy intensities relative to 2005 DO 5 IS=1,NS DO 5 IF=1,NF DO 5 IR=1,NR ENEW1(IS,IF,IR)=ENEW1A(IS,IF,IR)/EINT(IS,IF,IR) 5 EREN1(IS,IF,IR)=EREN1A(IS,IF,IR)/EINT(IS,IF,IR) C Read in file with future floor areas OPEN(unit=15,file='FloorAreaLow.data') c OPEN(unit=15,file='FloorAreaHigh.data') DO 10 I=1,3 10 READ(15,*) C The first data line is for 2005. The next 19 are at 5-yr intervals to 2100 DO 12 IY=0,19 12 READ(15,*) IYEAR,((AREAD(IS,IR,IY),IR=1,NR),IS=1,NS) C Interpolate floor areas at 5-yr intervals to floor areas every year DF=1.0/FLOAT(NY1) YEAR(0)=2005 DO 13 IS=1,NS DO 13 IR=1,NR AREAF(IS,IR,0)=AREAD(IS,IR,0) DO 14 I=1,19 DO 14 J=1,5 IY=J+5*(I-1) YEAR(IY)=YEAR(IY-1)+1 AREAF(IS,IR,IY)= * (AREAD(IS,IR,I-1)*FLOAT(5-J)+AREAD(IS,IR,I)*FLOAT(J))/5.0 C Net addition of floor area or net reduction of floor area each year ANEW(IS,IR,IY)=MAX(AREAF(IS,IR,IY)-AREAF(IS,IR,IY-1),0.0) 14 ARED(IS,IR,IY)=MAX(AREAF(IS,IR,IY-1)-AREAF(IS,IR,IY),0.0) C Divide the existing building stock into NY1 cohorts that will be either C replaced, renovated or skipped over each year from 2006 to 2050 inclusive DO 13 IY=1,NY1 AREPL(IS,IR,IY)=FREPL(IS,IR)*DF*AREAF(IS,IR,0) ARENO(IS,IR,IY)=FRENO(IS,IR)*DF*AREAF(IS,IR,0) AREMN(IS,IR,IY)=(1.0-FREPL(IS,IR)-FRENO(IS,IR))*DF*AREAF(IS,IR,0) C If floor area declines, ANEW is automatically zero but we need to reduce C the areas subject to renovation or replacement IF(ARED(IS,IR,IY).GT.0.0) THEN ARENO(IS,IR,IY)=ARENO(IS,IR,IY)-ARED(IS,IR,IY) IF(ARENO(IS,IR,IY).LT.0.0) THEN AREPL(IS,IR,IY)=AREPL(IS,IR,IY)+ARENO(IS,IR,IY) ARENO(IS,IR,IY)=0.0 IF(AREPL(IS,IR,IY).LT.0.0) WRITE(6,*) iy,is,ir,' arepl < 0' END IF END IF ATYEAR(IS,IR,IY)=AREPL(IS,IR,IY)+ARENO(IS,IR,IY) * + AREMN(IS,IR,IY)+ANEW(IS,IR,IY) 13 CONTINUE C Set the initial energy intensity of buildings that are renovated, or the energy C Divide existing stock into 45 cohorts and assign an initial energy intensity C distribution while still preserving the same mean intensity WRITE(6,140) WRITE(6,141) WRITE(6,142) DO 18 IR=1,NR 18 WRITE(6,143) REGION(IR),((EINT(IS,IF,IR),ENEW1A(IS,IF,IR), * EREN1A(IS,IF,IR),IF=1,NF),IS=1,NS) WRITE(6,*) DO 20 IS=1,NS DO 20 IF=1,NF WRITE(6,100) SECTOR(IS),FUEL(IF) WRITE(6,102) DO 20 IR=1,NR WRITE(6,104) REGION(IR),IYR0(IS,IF,IR),IYR1(IS,IF,IR), * ENEW0(IS,IF,IR),ENEW1(IS,IF,IR),EREN0(IS,IF,IR),EREN1(IS,IF,IR) DEVN=1.0-ENEW0(IS,IF,IR) DO 20 IY=1,NY1 C Deviation of the most-energy intensive and least energy intensive buildings C from the stock average. DEVN is chosen such that the last new buildings C built match the energy intensity of the least energy-intensive buildings of C the existing stock EISTRT(IS,IF,IR,IY)= EINT(IS,IF,IR)* * (1.0+2.0*DEVN*DF*(0.5*FLOAT(NY1)+0.5-FLOAT(IY))) 20 ETSTRT(IS,IF,IR,IY)=EISTRT(IS,IF,IR,IY)*AREAF(IS,IR,0)*DF C Verify that average energy initial energy intensity is unaltered (ie., that E0=1.0) c DO 25 IS=1,NS c DO 25 IF=1,NF c DO 25 IR=1,1 c E0=0.0 c DO 26 IY=1,NY1 c 26 E0=E0+EISTRT(IS,IF,IR,IY) c 25 WRITE(6,*) 'INITIAL ENERGY USE',E0*DF/EINT(IS,IF,IR) C Energy intensity of new and replacement buildings, first round C It is assumed that buildings are renovated or replaced in order of decreasing C energy intensity (this would happen if older buildings the most energy intensive C and buildings cohorts are renovated in order of decreasing age) DO 30 IS=1,NS DO 30 IF=1,NF DO 30 IR=1,NR N1=IYR0(IS,IF,IR)-ISTART N2=IYR1(IS,IF,IR)-ISTART DO 32 IY=1,N1 32 ENEW(IS,IF,IR,IY)=ENEW0(IS,IF,IR)*EINT(IS,IF,IR) DO 33 IY=N1+1,N2 33 ENEW(IS,IF,IR,IY)=EINT(IS,IF,IR)*(ENEW1(IS,IF,IR)*FLOAT(IY-N1) * +ENEW0(IS,IF,IR)*FLOAT(N2-IY))/FLOAT(N2-N1+1) DO 34 IY=N2+1,NY1 34 ENEW(IS,IF,IR,IY)=ENEW(IS,IF,IR,IY-1) DO 35 IY=1,NY1 35 EREPL(IS,IF,IR,IY)=ENEW(IS,IF,IR,IY) C Energy intensity of renovated buildings, first round DO 36 IY=1,N1 36 EREN(IS,IF,IR,IY)=EREN0(IS,IF,IR)*EINT(IS,IF,IR) DO 37 IY=N1+1,N2 37 EREN(IS,IF,IR,IY)=EINT(IS,IF,IR)*(EREN1(IS,IF,IR)*FLOAT(IY-N1) * +EREN0(IS,IF,IR)*FLOAT(N2-IY))/FLOAT(N2-N1+1) DO 38 IY=N2+1,NY1 38 EREN(IS,IF,IR,IY)=EREN(IS,IF,IR,IY-1) C Total energy consumption (TWh/yr) of buildings that are added, renovated, C replaced, or eligible for renovation or replacement, in year IY DO 30 IY=1,NY1 ETYEAR(IS,IF,IR,IY)= AREPL(IS,IR,IY)*EREPL(IS,IF,IR,IY) * +ARENO(IS,IR,IY)*EREN(IS,IF,IR,IY) * +AREMN(IS,IR,IY)*EISTRT(IS,IF,IR,IY) * +ANEW(IS,IR,IY)*ENEW(IS,IF,IR,IY) C Energy intensity (kWh/m2/yr) of the above buildings 30 EIYEAR(IS,IF,IR,IY)=ETYEAR(IS,IF,IR,IY)/ATYEAR(IS,IR,IY) C To compute the energy use in any given year, sum over all of the cohorts with C partial replacement or renovation, and all cohorts not yet touched DO 40 IS=1,NS DO 40 IF=1,NF DO 40 IR=1,NR ETOTAL(IS,IF,IR,0)=EINT(IS,IF,IR)*AREAF(IS,IR,0) DO 40 IY=1,NY1 ETOTAL(IS,IF,IR,IY)=0.0 DO 41 J=1,IY 41 ETOTAL(IS,IF,IR,IY)=ETOTAL(IS,IF,IR,IY)+ETYEAR(IS,IF,IR,J) DO 40 J=IY+1,NY1 40 ETOTAL(IS,IF,IR,IY)=ETOTAL(IS,IF,IR,IY)+ETSTRT(IS,IF,IR,J) C ROUND TWO C Set floor areas (billions m2) DO 50 IS=1,NS DO 50 IR=1,NR DO 50 IY=NY1+1,2*NY1 AREPL(IS,IR,IY)=FREPL(IS,IR)*ATYEAR(IS,IR,IY-NY1) ARENO(IS,IR,IY)=FRENO(IS,IR)*ATYEAR(IS,IR,IY-NY1) AREMN(IS,IR,IY)= * (1.0-FREPL(IS,IR)-FRENO(IS,IR))*ATYEAR(IS,IR,IY-NY1) C If floor area declines, ANEW is automatically zero but we need to reduce C the areas subject to renovation or replacement IF(ARED(IS,IR,IY).GT.0.0) THEN ARENO(IS,IR,IY)=ARENO(IS,IR,IY)-ARED(IS,IR,IY) IF(ARENO(IS,IR,IY).LT.0.0) THEN AREPL(IS,IR,IY)=AREPL(IS,IR,IY)+ARENO(IS,IR,IY) ARENO(IS,IR,IY)=0.0 IF(AREPL(IS,IR,IY).LT.0.0) WRITE(6,*) is,ir,iy,' arepl < 0' END IF END IF ATYEAR(IS,IR,IY)=AREPL(IS,IR,IY)+ARENO(IS,IR,IY) * + AREMN(IS,IR,IY)+ANEW(IS,IR,IY) 50 CONTINUE C Total energy consumption (kWh/yr) of buildings added or altered in year IY DO 52 IS=1,NS DO 52 IF=1,NF DO 52 IR=1,NR DO 52 IY=NY1+1,2*NY1 EISTRT(IS,IF,IR,IY)=EIYEAR(IS,IF,IR,IY-NY1) ETSTRT(IS,IF,IR,IY)=ETYEAR(IS,IF,IR,IY-NY1) C Energy intensities after second round of renovations EREPL(IS,IF,IR,IY)=EREPL(IS,IF,IR,IY-1)*(1.0-RATE) EREN(IS,IF,IR,IY)=EREN(IS,IF,IR,IY-1)*(1.0-RATE) ENEW(IS,IF,IR,IY)=ENEW(IS,IF,IR,IY-1)*(1.0-RATE) ETYEAR(IS,IF,IR,IY)= AREPL(IS,IR,IY)*EREPL(IS,IF,IR,IY) * +ARENO(IS,IR,IY)*EREN(IS,IF,IR,IY) * +AREMN(IS,IR,IY)*EISTRT(IS,IF,IR,IY) * +ANEW(IS,IR,IY)*ENEW(IS,IF,IR,IY) 52 EIYEAR(IS,IF,IR,IY)=ETYEAR(IS,IF,IR,IY)/ATYEAR(IS,IR,IY) C Compute total stock energy use for each year in round two DO 54 IS=1,NS DO 54 IF=1,NF DO 54 IR=1,NR DO 55 IY=NY1+1,2*NY1 ETOTAL(IS,IF,IR,IY)=0.0 DO 56 J=NY1+1,IY 56 ETOTAL(IS,IF,IR,IY)=ETOTAL(IS,IF,IR,IY)+ETYEAR(IS,IF,IR,J) DO 55 J=IY+1,2*NY1 55 ETOTAL(IS,IF,IR,IY)=ETOTAL(IS,IF,IR,IY)+ETSTRT(IS,IF,IR,J) DO 54 IY=2*NY1+1,NY 54 ETOTAL(IS,IF,IR,IY)=ETOTAL(IS,IF,IR,2*NY1) C Compute sector, fuel, and global sums DO 60 IY=0,NY DO 61 IR=1,NR AREAF(NS1,IR,IY)=AREAF(1,IR,IY)+AREAF(2,IR,IY) DO 62 IS=1,NS 62 ETOTAL(IS,NF1,IR,IY)=ETOTAL(IS,1,IR,IY)+ETOTAL(IS,2,IR,IY) DO 61 IF=1,NF1 61 ETOTAL(NS1,IF,IR,IY)=ETOTAL(1,IF,IR,IY)+ETOTAL(2,IF,IR,IY) DO 60 IS=1,NS1 AREAF(IS,NR1,IY)=0.0 AREAF(IS,NR2,IY)=0.0 AREAF(IS,NR3,IY)=0.0 DO 63 IR=1,NR 63 AREAF(IS,NR+1,IY)=AREAF(IS,NR1,IY)+AREAF(IS,IR,IY) C Compute separate sums for OECD/EEur/FSU and other countries DO 64 IR=1,5 64 AREAF(IS,NR2,IY)=AREAF(IS,NR2,IY)+AREAF(IS,IR,IY) DO 65 IR=6,NR 65 AREAF(IS,NR3,IY)=AREAF(IS,NR3,IY)+AREAF(IS,IR,IY) DO 60 IF=1,NF1 ETOTAL(IS,IF,NR1,IY)=0.0 ETOTAL(IS,IF,NR2,IY)=0.0 ETOTAL(IS,IF,NR3,IY)=0.0 DO 66 IR=1,NR 66 ETOTAL(IS,IF,NR1,IY)=ETOTAL(IS,IF,NR1,IY)+ETOTAL(IS,IF,IR,IY) DO 67 IR=1,5 67 ETOTAL(IS,IF,NR2,IY)=ETOTAL(IS,IF,NR2,IY)+ETOTAL(IS,IF,IR,IY) DO 60 IR=6,NR 60 ETOTAL(IS,IF,NR3,IY)=ETOTAL(IS,IF,NR3,IY)+ETOTAL(IS,IF,IR,IY) C IR=NR1 C DO 70 IR=1,NR3 WRITE(6,120) REGION(IR) WRITE(6,121) (SECTOR(IS),(SECTOR(IS),IF=1,NF),IS=1,NS1), * ((SECTOR(IS),IF=1,NF),IS=1,NS) WRITE(6,122) (AREAN,(FUEL(IF),IF=1,NF),IS=1,NS1), * ((FUEL(IF),IF=1,NF),IS=1,NS) WRITE(6,123) (UNITS(1),(UNITS(2),IF=1,NF),IS=1,NS1), * ((UNITS(3),IF=1,NF),IS=1,NS) DO 70 IY=0,NY,5 70 WRITE(6,124) YEAR(IY), * (AREAF(IS,IR,IY),(ETOTAL(IS,IF,IR,IY),IF=1,NF),IS=1,NS1), * ((ETOTAL(IS,IF,IR,IY)/AREAF(IS,IR,IY),IF=1,NF),IS=1,NS) C STOP C Detailed outputs for analysis purposes DO 80 IS=1,NS DO 80 IF=1,2 DO 80 IR=1,3 WRITE(7,130) REGION(IR),SECTOR(IS),FUEL(IF) WRITE(7,131) DO 81 IY=1,NY1 81 WRITE(7,132) IY,AREPL(IS,IR,IY),ARENO(IS,IR,IY),AREMN(IS,IR,IY), * ANEW(IS,IR,IY),ATYEAR(IS,IR,IY),EISTRT(IS,IF,IR,IY), * ENEW(IS,IF,IR,IY),EREN(IS,IF,IR,IY),ETYEAR(IS,IF,IR,IY), * ETSTRT(IS,IF,IR,IY),EIYEAR(IS,IF,IR,IY),ETOTAL(IS,IF,IR,IY) WRITE(7,131) DO 80 IY=NY1+1,NY 80 WRITE(7,132) IY,AREPL(IS,IR,IY),ARENO(IS,IR,IY),AREMN(IS,IR,IY), * ANEW(IS,IR,IY),ATYEAR(IS,IR,IY),EISTRT(IS,IF,IR,IY), * ENEW(IS,IF,IR,IY),EREN(IS,IF,IR,IY),ETYEAR(IS,IF,IR,IY), * ETSTRT(IS,IF,IR,IY),EIYEAR(IS,IF,IR,IY),ETOTAL(IS,IF,IR,IY) WRITE(8,150) SECTOR(1) WRITE(8,150) (((REGION(IR),I=1,3),IF=1,2),IR=1,4,3) WRITE(8,150) (((FUEL(IF),I=1,3),IF=1,2),IR=1,4,3) WRITE(8,150) (((VARIAB(I),I=1,3),IF=1,2),IR=1,4,3) DO 90 IY=1,NY 90 WRITE(8,152) IY,((EISTRT(1,IF,IR,IY),ENEW(1,IF,IR,IY), * EREN(1,IF,IR,IY),IF=1,2),IR=1,4,3) WRITE(8,150) (((REGION(IR),I=1,3),IF=1,2),IR=1,4,3) WRITE(8,150) (((FUEL(IF),I=1,3),IF=1,2),IR=1,4,3) WRITE(8,150) (((VARIAB(I),I=1,3),IF=1,2),IR=1,4,3) C 100 FORMAT(' ',2A10) 102 FORMAT(' Region IYR0 IYR1 ENEW0 ENEW1 EREN0 EREN1') 104 FORMAT(' ',A10,2I6,4F8.2) 120 FORMAT(/,' ',A10,' region') C -------!---------!---------!---------!---------!---------!---------!-- 121 FORMAT(' - ',13A10) 122 FORMAT(' - ',13A10) 123 FORMAT(' - ',13A10) 124 FORMAT(I6,13F10.2) 130 FORMAT(3A10) 131 FORMAT(' AREPL AREN AREMN ANEW ATYEAR EISTRT ENEW * EREP ETYEAR ETSTRT EIYEAR ETOTAL') 132 FORMAT(I4,5F7.3,5F9.3,2F10.3) 140 FORMAT(12X,' - Res Res Res Res Res Res Com *m Comm Comm Comm Comm Comm') 141 FORMAT(12X,' - Elec Elec Elec Fuels Fuels Fuels Ele *c Elec Elec Fuels Fuels Fuels') 142 FORMAT(4X,' Region Start Finish Finish Start Finish Finish * Start Finish Finish Start Finish Finish') 143 FORMAT(' ',A10,12F8.1) 150 FORMAT(7X,12(1X,A10)) 152 FORMAT(' ',I6,12F11.3) C STOP END